Quantum tunneling as a classical anomaly 



Carl M Bender* and Daniel W Hook** 

*Department of Physics, Washington University, St. Louis, MO 63130, USA 
email: cmb@wustl.edu 

t Theoretical Physics, Imperial College, London SW7 2AZ, UK 
email: d.hook@imperial.ac.uk 

Abstract. Classical mechanics is a singular theory in that real-energy classical 
particles can never enter classically forbidden regions. However, if one regulates 
classical mechanics by allowing the energy E of a particle to be complex, the particle 
exhibits quantum-like behavior: Complex-energy classical particles can travel between 
classically allowed regions separated by potential barriers. When ImE — » 0, the 
classical tunneling probabilities persist. Hence, one can interpret quantum tunneling as 
an anomaly. A numerical comparison of complex classical tunneling probabilities with 
quantum tunneling probabilities leads to the conjecture that as Re E increases, complex 
classical tunneling probabilities approach the corresponding quantum probabilities. 
Thus, this work attempts to generalize the Bohr correspondence principle from 
classically allowed to classically forbidden regions. 
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Classical mechanics is a singular limit of quantum mechanics. (The limit e — > 
is singular if an abrupt change occurs at e = 0. For example, lim^o e x 5 + x — 1 
is singular because four of its five roots abruptly disappear at e = 0.) The classical 
limit h — >■ of the time-independent Schrodinger equation h 2 ip"(x) = [V(x) — E]ip(x) 
is singular because at H — it is no longer possible to impose initial or boundary 
conditions on the wave function ip{x). Moreover, while quantum particles are able to 
enter classically forbidden regions, these particles abruptly lose this ability at H — 0. 
Thus, the phenomenon of tunneling seems to be entirely quantum mechanical, and 
one may not ask such classical questions as, Which path does the particle follow while 
tunneling? 

It may be possible to recover some features of a theory that were abruptly lost in a 
singular limit, like the ability of a particle to enter a classically forbidden region. To do 
so, one introduces a regulator, which is then removed in a careful limiting process. The 
features that remain after the regulator is removed are referred to as anomalies. (The 
axial anomaly in quantum field theory can be obtained by using dimensional regulation.) 

Particle trajectories in conventional classical mechanics are real functions of time, 
but recent studies of the complex solutions to the classical equations of motion 
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(Hamilton's equations) have shown that real-energy classical particles may leave the real 
axis and travel through complexified coordinate space [TjISlEllIlElEJlTllHllH], [10l [HI [12] . 
While these complex classical trajectories may pass through classically forbidden regions 
on the real axis, complex classical mechanics is still a singular theory because there is 
no tunneling; that is, no complex path runs from one classically allowed region on the 
real axis to another. 

Complex classical mechanics may be regulated by taking the energy of a particle 
to be complex. Surprisingly, a complex-energy classical particle exhibits qualitative 
features normally associated with a quantum particle: Such a classical particle can 
exhibit tunneling-like behavior in which it travels from one classically allowed region to 
another classically allowed region even though these two regions are separated on the real 
axis by a classically forbidden region [13, 14J. Furthermore, a complex-energy classical 
particle in a periodic potential exhibits quantum-like behavior; there are sharply defined 
energy bands separated by gaps. In these energy bands the classical particle exhibits 
a kind of resonant tunneling [13] . The time-energy uncertainty principle supports our 
choice of regulator. This uncertainty principle implies that a precise measurement of 
the energy of a particle in a finite time interval is impossible; some uncertainty AE is 
associated with such a measurement. In this paper we allow AE to be complex. 

This paper explores the connection between standard quantum mechanics and 
complex-energy classical mechanics at a quantitative level. We show that as the regulator 
is removed (ImE — > 0), the tunneling behavior of complex-energy classical mechanics 
persists. Thus, we recover tunneling as a kind of anomaly. 

To compare the behavior of quantum particles and complex classical particles, we 
use quartic and sextic asymmetric double-well potentials and compute numerically the 
relative probabilities of finding these particles in each well. We find that as the number 
of nodes in the quantum-mechanical eigenfunction increases, the quantum and classical 
probabilities approach one another. Thus, we demonstrate that the Bohr correspondence 
principle, which has recently been generalized to the complex domain [151 [16], actually 
applies to tunneling phenomena. 

The quartic double-well potential 



H = p 2 + V^(x), where we have set U = 1, are Eq = —18.0182 (below the bottom 
of the right well), E 1 = -7.1879, E 2 = -6.8595, E 3 = 1.6806, and E A = 2.8845 (above 
the bottom of the right well and below the peak of the barrier, and E$ = 8.3312 (above 
the barrier). 

We plot the sextic potential 



h = 1 is Eq = —23/2. (This exact value of Eq provides a benchmark confirming that our 
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Figure 1. Upper panel: Quartic asymmetric double-well potential V^\x) in |l| 
and the first six quantum energy levels. There are minima at x = —1.2499, where 
= -24.0384, and at x = 1.9165, where = -12.5501. The peak of the 

barrier is at x — 0.4884, where — 4.1144. The ground-state energy Eq lies below 
the bottom of the right potential well. The next four energy levels lie between the 
bottom of the right potential well and the top of the barrier. The sixth energy level E$ 
lies above the barrier. Lower panel: Sextic asymmetric double-well potential (x) 
in The minima of the wells are at x = -1.7083, where V {6) = -20.7710, and at 
x = 1.8215, where = -13.9373. The peak of the barrier is located at x = -0.5184, 
where = 4.2731. 
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Figure 2. Quantum probability density for a particle in the nth eigenstate in the 
quartic potential V^(x) in (jlj). The probability P^^ n in (J3J of finding a particle 
in the nth eigenstate to the right of the barrier is determined by integrating the nth 
probability density. 



numerical calculations are accurate to better than 13 decimal places.) The first three 
excited quantum states have energies E\ = —9.9690, E 2 = —3.9819, and E% = 1.8095. 

The eigenfunction ip n (x) corresponding to the energy E n has n nodes. In Fig. [2] we 
plot the probability density |?/>n(£)| 2 for a particle in V^'(x) for n = 0, 1, 5. We 
integrate \il) n {x)\ 2 to determine the probability P^Tn °f finding the particle to the right 
of the barrier 
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(3) 
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For the sextic potential V^ 6 \x) in ^ we find that for the first four quantum states the 
probabilities of finding the particle to the right of the top of the barrier are 

P%$i> = 0.0391%, = 99.9986%, 

p right n 2 = 99.8651%, P r ^ 3 = 78.7223%. (4) 
We now describe the behavior of a classical particle in the quartic potential Q for 
the case in which x is complex. Such a particle may have either real or complex energy, 
and we begin by considering the complex trajectories of a real-energy classical particle. 
Let us take the energy of a classical particle in to be E4 = 2.8845. There are two 
real turning points at x — —1.9428 and 0.2234, which bound the classically allowed 
region to the left of the barrier and two turning points at x = 0.7596 and 2.4998, which 
bound the classically allowed region to the right of the barrier. If the initial position 
of the particle is real and in the classically allowed region in the right (or left) well, 
the trajectory of the particle oscillates on the real axis between the right (or left) pair 
of turning points. However, if the particle is initially in a classically forbidden region 
on the real axis, the particle leaves the real axis and travels in a closed periodic orbit 
in the complex-x plane [see Fig. [3] (upper panel)]. The orbits in the complex-x plane 
enclose the classically allowed regions on the real axis and never cross the vertical line 
Rex = 0.4884, which passes through the peak of the barrier. The periods T of all closed 
orbits, both on the left and right side of Fig. |3j are the same and are given by 

fright fa 

= y, lcft v/p-v^h*)' 

where x\ c tt and a^ght are the left and right turning points in either well and E is the 
(real) classical energy of the particle. For the case in which the energy E = E^ = 2.8845, 
the period T = 0.8464. The particle in the upper panel of Fig. [3] cannot travel from one 
classically allowed region to the other classically allowed region, so there is no tunneling 
effect when the classical energy is exactly real. 

We now explain how deterministic classical systems can produce results that are 
numerically comparable to the quantum results in ^ and Q. If the energy E of 
the classical particle is complex, the classical trajectory of such a particle is not in 
general closed [13] . In Fig. [4] we plot the complex path of a classical particle of energy 
E = 2.8845 + 0.25i for the time interval from t = to t = T = 0.8464. Observe that 
while the particle executes a loop of nearly 360° around the two turning points, the 
trajectory of the particle is no longer closed. If this particle is initially in a classically 
allowed region, it spirals outward around the pair of turning points that bound the 
region. However, the particle does not drift off to infinity. Rather, it crosses the vertical 
line that passes through the top of the barrier and spirals into the other well. This 
behavior is shown in Fig. [3] (lower panel), where a particle of energy E = 2.8845 + 0.5i 
alternately visits both potential wells. This classical particle executes deterministic 
tunneling from well to well. 

There is an unexpected subtlety in the behavior of complex-energy classical particles 
that should not be overlooked. While the spiral motion shown in the lower panel in Fig. [3] 
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Figure 3. Upper panel: Periodic trajectories in the complex- x plane of a classical 
particle of real energy E^ = 2.8845 in the quartic potential 0. The dots indicate 
the four turning points. There is no tunneling-like behavior; no trajectory runs from 
one classically allowed region to the other. Lower panel: Classical trajectory in the 
complex- a; plane of a particle of complex energy E = 2. 8845 + 0. 5z in the same potential 
Unlike the closed trajectories in the upper panel, this classical trajectory is open. 
The particle begins at x = 1 and traces an outward anticlockwise spiral around the 
right pair of turning points. It then crosses to the left of the peak of the barrier and 
spirals inward and clockwise around the left pair of turning points. Next, it crosses the 
real axis between the left turning points and spirals outward and anticlockwise around 
the left pair of turning points. The time for the displayed trajectory is t = 8. For 
t > 8, the particle continues to spiral inward and outward as it oscillates from well to 
well. The particle exhibits this deterministic tunneling-like behavior for all time, but 
the trajectory never crosses itself. This trajectory behaves as if it were controlled by 
a pair of strange attractors. 




-2 L 



Figure 4. Complex path of a particle of complex energy E — 2.8845 + 0.25i. The 
time of this path is taken to be t = T = 0.8464, which is the same as the period of a 
classical particle of real energy E — 2.8845. Note that the particle in this figure makes 
an approximately full angular revolution around the two two turning points, which are 
denoted by dots. However, the path is not closed. 



is typical of complex-energy classical particles, it has recently been discovered that for 
a dense set of measure zero of complex classical energies the classical paths are periodic 
[2]. The complex energies for which the classical motion remains closed and periodic 
lie on an infinite number of nearly straight lines that eminate from the origin in the 
complex-energy plane. A schematic representation of these lines is shown in Fig. [5] 

The key advance in this paper is the numerical observation that the probabilities 
associated with the tunneling of a classical particle having complex energy persist as 
the regulator ImAi? tends to zero even though there is no classical tunneling when 
ImE = 0. The upper panel in Fig. [3] shows that a classical particle having real energy 
does not exhibit tunneling behavior. Yet, a classical particle having ImE ^ typically 
does exhibit tunneling behavior, and the corresponding probabilities can be compared 
in the limit as ImE — > with the quantum tunneling probabilities in ^ and Q. 

To demonstrate the persistence of classical tunneling in the limit ImE — > 0, we 
compute the classical trajectory for long times [17] and determine the fraction of time 
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Figure 5. Schematic drawing of the special complex energies for which the particle 
motion in the quartic potential ([!]) is periodic for any choice of initial condition. These 
energies lie on asymptotically straight dashed curves that approach the origin in the 
complex-£? plane. The energies form a set of measure zero; there are an infinite 
number of such curves and these curves are dense in the complex-_E plane. The 
regulation procedure used in this paper establishes a connection between complex 
classical mechanics and quantum mechanics by taking a sequence of complex classical 
energies in which the imaginary part tends to 0. This sequence is indicated by dots 
on the solid vertical line. Although this line intersects the dashed lines infinitely many 
times, the numerical procedure remains unaffected. 

that the classical particle spends to the left and to the right of the line Re x = 0.4884 (the 
location of the peak of the potential barrier). As ImE gets smaller, the fraction of time 
spent to the right (and to the left) of the barrier approaches a constant]!] Take the real 
part of the classical energy to be the fourth excited quantum energy Re E = E± = 2.8845 
and take ImE = 2~ k . For k — 0, 1, 2, 3, 4, 5 the chance of finding the classical particle 
to the right of the barrier is 50.4%, 53.9%, 55.0%, 53.1%, 52.8%, 53.3%. This sequence 
approaches a limiting anomalous value of about 53%, which exceeds 50% in accordance 
with the heuristic semiclassical argument in Ref. [18]. A classical version of the time- 

| Figure [4] shows that each angular revolution in the complex-energy spiral takes approximately the 
same time regardless of whether the spiral encircles the left or the right pair of turning points. The 
left-right asymmetry in the well gives rise to different left and right winding numbers in the spiral path. 
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energy uncertainty principle applies here: As the imaginary part of the energy gets 
smaller, it takes more time for the particle to oscillate between the wells (13], and 
thus this extrapolation procedure for determining the classical tunneling probabilities 
requires more computer time. 

Classical tunneling persists in the limit as ImE — > even though there is no 
classical tunneling at ImE = 0; Fourier series exhibit strongly analogous behavior. 
Take the function f(x) = 1 on [0, ir]. Although /(0) ^ and f(Tr) ^ 0, we can 
represent f(x) as the Fourier sine series | ^2^=0 2fc+l sin[(2£; + l)x}. The partial sum 
Sk(x) = J2k=o ak sm (kx) converges to f{x) as K — )■ oo on the open interval (0, n) 
because f(x) is continuous. All terms in Sk(x) vanish at and n, but we can still recover 
the nonzero value of /(0) and f(ir) from Sk(x). The extrapolation procedure used 
above for classical tunneling probability can be applied to Fourier series to determine 
/(0). We take twice as many terms in Sk{x) as x is halved and thereby circumvent 
the problem of nonuniform convergence (the Gibbs phenomenon): We evaluate Sk(x) 
at K — 100 x 2 k and x = 2~ k for k = 0, ...,5. The numerical values of Sk{x) 
are S 100 = 0.997776, S 200 = 0.996 704, S 400 = 0.997293, S mo = 0.997818, S imo = 
0.998 128, S 3200 = 0.998 292. We infer that /(0) = 1. 

In analogy with the extrapolation scheme used for Fourier series, we determine 
the classical probability of finding the particle to the right of the peak of the barrier 
for a sequence of energies in which the real part of the energy is held fixed and the 
imaginary part of the energy tends to zero. For a small but fixed imaginary classical 
energy, the classical probability approaches the quantum probability as the real part 
of the energy increases. For example, when ImE' = 1/4 and KeE = E n , where E n 
is the nth quantum eigenenergy for the quartic potential V^(x) in Q, we obtain the 
following classical probabilities for finding the classical particle in the right well: 

pclass KK 4% pclass _ rr Q(W 

-'right, 1 — OO.^/O, "'rights — OO.U/0, 

P&W = 54-3%, P r f g a h s t s ,4 = 55.0%. (5) 

The quantum probabilities P^^n m oscillate about these classical probabilities and 
are in good agreement when n = 3 and n — 4. For deeper double-well potentials, the 
classical probabilities continue to approach the quantum probabilities as n increases. 

For the sextic potential V^ 6 \x) in ([2J there are six turning points. When the 
classical energy is real and between the minimum of the lower well and the top of the 
barrier, the turning points group into three pairs, one pair on the real axis to the left 
of the barrier, a second pair on the real axis to the right of the barrier, and a complex- 
conjugate pair that is associated with the complex extension of the barrier. If we take 
the classical energy to be real, say E n , we cannot observe tunneling because the classical 
orbits are closed and periodic and just encircle the pairs of turning points, as shown in 
the upper panel of Fig. |6j 

When the classical energy E is complex, the classical path is no longer closed. A 
particle trajectory beginning at x — 1 in the right well spirals outward and eventually 
circles around the barrier turning points. The trajectory does not always penetrate 
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Figure 6. Complex classical trajectories for a particle subject to the sextic potential 
(|]). The energy of the particle in the upper panel is real and three closed 
nonintersecting trajectories for this particle are shown. The energy of the particle 
in the lower panel is complex and a single open trajectory is plotted. 



to the left well; sometimes the particle is ejected and falls back into the right well. 
However, after sufficiently many tunneling attempts the particle spirals inward around 
the left pair of turning points, as shown in the lower panel of Fig. [6] Taking Im E = 1/16 
and Re E = E n (n = 0, 1, 2, 3), we find that as the real part of the energy increases, the 
classical tunneling probabilities i^ghtn listed below approach the quantum probabilities 
^£ in 

pclass _ qi 70/ pclass _ 07 i 07 

-^right.O - - ' /0, -^right.l —O(.l/0, 

P&* = 32-4%, P r f g a h s t s , 3 = 72.3%. (6) 
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The agreement between Mj) and ([6]) is better than that for the quartic case, possibly 
because the mixture of successful and failed tunneling attempts through the central 
barrier region leads to a more accurate analog of quantum tunneling. We conclude 
that complex classical mechanics provides a good approximation to quantum tunneling 
for the higher-energy states, as one would expect of a generalized Bohr correspondence 
principle. 

Of course, the agreements between (|3| and (]5]), and (|4]) and (|6| are not very precise. 
To obtain better agreement, it is necessary to take deeper asymmetric double wells 
having more quantum eigenenergies below the top of the central barrier and above 
the bottom of the shallower well. This can also be achieved by taking h to have a 
smaller numerical value. (The coefficient of —d 2 /dx 2 in the Schrodinger equation is 
H 2 .) Until now, we have taken h = 1 but we have repeated our numerical work for 
the quartic potential ([I]) with h = 1/2. Now, the highest eigenvalue below the top of 
the barrier is E 8 = 3.7239. The corresponding probability density iV'sf^)! 2 is shown in 
Fig. [7j For this eigenf unction, the quantum probability of finding the particle to the 
right of the peak of the barrier is 60.29% and the corresponding classical probability 
using the complex energy E = E 8 + i/4 is 55.2%. (These numbers compare favorably 
with -Prigh t nt 4 = 40.7689% and ^^,4 = 55.0% for the highest eigenvalue in the double 
well with h — 1.) Doing further calculations for successively smaller values of h would 
require a considerable increase in computer time. The purpose of this paper is only 
to demonstrate, in principle, our conjectured connection between complex classical 
and quantum tunneling, and not to carry out an extensive and detailed numerical 
investigation. 

We thank the U.S. Department of Energy for financial support and Mathematica 
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Figure 7. Quantum probability density for a particle in the i/js( x ) eigenstate in the 
quartic potential with ft = 1/2. The energy of the particle is Eg = 3.7239. The 
probability density has eight nodes. Note that the particle is most likely to be found 
inside or near the classically forbidden region associated with the barrier; this shows 
that the particle does not have a high enough quantum number for its behavior to 
resemble that of a classical particle. 
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